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1.0  INTRODUCTION 

The  objective  of  the  proposed  research  is  to  develop  predictive  models  characterizing  the 
response  due  to  deformation,  diffusion  and  reaction  in  service  environments  for  PMR-15  resin 
composites.  The  principal  investigator  has  developed  a  theoretical  model  from  a  fully 
thermodynamic  standpoint  to  describe  the  response  of  polymeric  composites  that  are  subject  to 
diffusion  and  chemical  reactions  due  to  hostile  environmental  conditions.  The  thermodynamic 
framework  is  very  general  and  has  been  used  to  describe  the  disparate  and  diverse  response  by  a 
large  class  of  materials:  viscoelastic  fluids  (Rajagopal  and  Srinivasa  (2000),  Krishnan  and 
Rajagopal  (2004)),  inelasticity  (Rajagopal  and  Srinivasa  (1998)),  crystallization  of  polymers 
(Rao  and  Rajagopal  (2001),  twinning  (Rajagopal  and  Srinivasa  (1997)),  Kannan  et  al.  (2002)), 
shape  memory  alloys  (Rajagopal  and  Srinivasa  (1999)),  single  crystal  super  alloys  (Prasad  et  al. 
(2005),  Prasad  et  al.  (2006)),  and  viscoelastic  solids  (Kannan  and  Rajagopal  (2004)).  There  are 
two  central  ideas  on  which  the  framework  is  built.  The  first  is  that  when  a  body  undergoes  an 
entropy  producing  process,  the  underlying  natural  configuration"  of  the  body  (see  Rajagopal 
(1995)  for  a  discussion  of  natural  configuration)  usually  changes  and  the  second,  the  constitutive 
relation  that  describes  the  response  of  a  body,  from  amongst  a  competing  class  of  possible 
response  functions  is  one  that  maximizes  the  rate  at  which  entropy  is  produced.  In  this  work, 
using  both  these  above  ideas,  a  model  is  developed  to  describe  the  response  of  a  viscoelastic 
polymeric  solid  due  to  its  deformation  when  a  fluid  is  diffusing  through  the  solid,  taking  into 
account  temperature  effects.  Such  a  model  can  be  used  to  describe,  in  addition  to  PMR-15  resin 
composite,  a  whole  host  of  viscoelastic  polymer  composites  due  to  hostile  environmental 
conditions,  namely  high  temperature  and  the  diffusion  of  moisture  through  the  body. 
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2.0  DEVELOPMENT  OF  A  VISCOELASTIC  SOLID  MODEL  WITH  DEGRADATION 
DUE  TO  TEMPERATURE 

In  this  section,  we  shall  first  develop  a  model  for  the  viscoelastic  behavior  of  PMR-15.  This 
model  also  accounts  for  degradation  due  to  high  temperatures. 

2.1  Preliminaries 


Let  kr  be  the  reference  configuration  and  nt  denote  the  current  configuration.  Then,  the  motion 
Xkr  is  the  one  -  one  mapping  that  assigns  to  each  point  X  €  kr,  a  point  x  e  nt,  at  a  time  t,  i.e., 

x  =  xKfi(X.t).  (i) 


Let  Kp(t)  denote  the  natural  configuration  (stress  -  free)  configuration  corresponding  to  Kt.  This 
natural  configuration  is  reached  instantaneously  from  the  current  configuration,  by  the  body, 
upon  removal  of  a  load.  The  deformation  gradient  F  is  defined  by 

F ,=  dX»H 

'  r>X  '  (2) 


Let  r  kp(o  be  the  gradient  of  the  mapping  from  kp( 0  and  let  G  be  the  gradient  of  mapping  from 
KR  to  kp(0  (see  Figure  1).  Then, 


F  =  F«*oG’ 


(3) 


or 


G 


F. 


We  shall  denote  the  left  Cauchy  -  Green  stretch  tensors 


Bg  :=  GGr, 


B 


p(t) 


:=  F. 


p(i)  Kp(t) ' 


We  shall  also  define  the  velocity  gradients 

L6  :=  GG-1,  L  =  FF  \ 


FS(()F«P(0’ 


with  their  symmetric  parts  by 

D,  =  -  (L,  +  Lf ) ,  i  =  p(t),  G  or  no  subscript. 

Also,  we  define 

lBi=tr(  B,),  I  In,  =  ^[(tr(Bl)f-tr(Bj)},  /  =  G,p(t), 

where  1 ' 1  •  ■  is  the  trace  operator  for  a  second  order  tensor.  Now, 

F  =  F^(0G  +  FS(0G 
=>  FF-1  =  FKp(t)GG-lF”p1(t)  +FKj>(t)G 

=>  L  =  Lp(t)  +  FKp(t)LcF~^t), 

where  t  • )  is  the  time  derivative  of  the  second  order  tensor.  In  addition, 


(4) 

(5) 

(6) 

(7) 

(8) 


(9) 
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where  ( •)  is  the  time  derivative  of  the  second  order  tensor.  In  addition, 


and  similarly 


Br(t)  =  F«mF 
=  Lp(t)Bp(t)  + 


Bg  =  LgBg  +  BgLg- 


Hence,  from  (9)  and  (10),  we  have 

B/>(t)  =  BBp(t)  +  Bp{t)LTp(t)  -  FKp(f)  (Lg  +  L£)  F^, 

and  so 


(10) 

(11) 

(12) 


v 


B  p{t)  2FKp(4)DgF 


Kp(t ) 


(13) 


v 

where  (  )  is  the  Oldroyd  derivative.  The  non  -  isothermal  reduced  energy  dissipation  equation 
(Kannan  and  Rajagopal  (2004))  is  given  by: 


T.D  —  pip  —  psO 


q/,  .grad  (9 
9 


=  P0(  :=  £  >  0, 


(14) 


where  T  is  the  Cauchy  stress,  '  ’  is  the  specific  Helmoltz  free  energy,  P  is  the  density  of 
viscoelastic  solid,  ^  is  the  temperature,  s  is  the  specific  entropy,  O/i  is  the  heat  flux,  C  is  the  rate 
of  entropy  production  and  £  is  the  rate  of  dissipation. 


2.2  Constitutive  assumptions  and  maximization  of  the  rate  of  dissipation 
2.2.1  General  results 


Now,  we  shall  assume  that  the  viscoelastic  solid  is  isotropic  and  incompressible  with  the  specific 
Helmoltz  potential  of  the  form 


Ip  -  1p(Bp(t),BG,9)  -  Tp(lBpW,IlBp(t)JBn<IlBn<9)- 


Jp(D’ 


Based  on  the  assumption  of  isotropic  elastic  response,  we  choose  Kp(f)  such  that 

F  —  V 

KP(t)  KP(t )’ 

V  •  •  F 

where  Kp<t>  is  the  right  stretch  tensor  in  the  polar  decomposition  of 

We  shall  also  assume  that  the  total  dissipation  can  be  split  additively  as  follows 

T.D  —  pip  —  p80  =  £m  >  0.  -fWSrad  ^  =  ^,  >  (K 

where  Pm,  Pc  are  the  rates  of  mechanical  dissipation  and  dissipation  due  to  conduction 
respectively.  Now,  we  constitutively  choose 

q/j  =  —k(0) grad#,  k(6)  >  0, 


(15) 

(16) 


(17) 


(18) 


where  k  is  the  thennal  conductivity;  so  that  (17(b))  is  automatically  satisfied. 
Now, 
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J,  _  |Y  ,  J  gj  ^  T  _  fo/i  1  • 

+  f  (#- + /b0J^)  i  - ,;';  h.  I  bg + §i 

\OIbg  dIlBGJ  9IIbg  96 


L  \  /  J 

and  using  ( 1 0),  ( 1 1 )  along  with  ( 1 6)  in  ( 1 9) 

_  2  |Y  <W'  r  9j>  \  _  _gg_  2 


.(D  -  D 


Also,  we  assume 


,  n\(  ,  r  \r  ^  n2l 

\9Ibg  9IIbg  J  9IIbg  J  d# 

he  mechanical  rate  of  dissipation  to  be  of  the  form 

Cm  =  Cm  ($  >  Bp(f )  >Dg). 


On  substituting  (20)  into  (17(a)) 

rr,  o  (  9$  ,  T  9rp  \T3  ,  0  9$  „2 

[  - fc  +  B",,) + 

,  n  /  94)  (  ,  d'tjj  n2 

^  [fe  /B""  ^2)  B"('> " 

o,  Y  ^  If  ^  ^  pal  n 

P  +  Ibg  91 1 Bg)  G  9IIBg  G\  G 

Ta.7.  1 


We  shall  choose 


9^  a 

-'’[w+y 

=  Cm(^,Bp(t),DG). 


and  define 


9ip 

S  ~  ~96' 


Tpm  p,I  +  2p  [(^  +  - 


where  Pi,  P2  are 


T0:=KI  +  2p[f^f +/Bo^)b0-t 

the  Lagrange  multipliers.  Using  (23)-(25)  in  (22; 

(T  -  Tp(t))  .D  +  (Tp(t)  -  To)  .D( 


(T  -  Tp(t))  .D  +  (Tp(, 
=  Cm(^Bp(t).DG). 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
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From  constraint  of  incompressibility,  we  have 

tr(  D)  =  tr(Dp{t))  =  tr{  DG)  =  0. 

(27) 

Since,  RHS  of  Eqn.(26)  does  not  depend  on  using  (27), 

T  =  pi  + 

(28) 

where  P  is  the  Lagrange  multiplier  due  to  incompressibility,  with 

(T p(t)  -  TG)  ,Dg  =  ^m(Bp((),DG), 

(29) 

which  can  be  re  -  written  as 

(T-TG).DG  =  £m(BpW,DG), 

(30) 

using  (27)  and  (28)  . 

Now,  following  Rajagopal  and  Srinivasa  (2000),  we  maximize  the  rate  of  dissipation  by 

varying  DG  for  fixed  For  this,  we  maximize  the  function 


4*  :=  U  +  Ai  [U  -  (T  -  Tg)  .Dg]  +  A2(I.Dg), 
where  Ai,  A2  are  the  Lagrange  multipliers.  By  setting,  ()<$>/ ODo  —  0,  we  get 

Ai  +  1  \  0£m 


A2 

t  =  tg  +  -^i  + 
M 


On  substituting  (32)  in  (30)  we  get 


Ai  +  1 

Ai 


Ai 

U 


dBG 


d£n 


and  so  (32)  with  (25)  becomes 

dip 


0DC 


D 


G 


T  =  2  p 


01  Da 


+  Ibg  ~ 


dip 


OIIbo 


bg- 


dip 

OIIBg 


B  % 


+ 


& 


9  in 


.  OH  a 


_ \  (Krn_ 

dg  )  £>dg 


+  AI. 


where  A  :=  ^  +  p2  is  the  Lagrange  multiplier  due  to  the  constraint  of  incompressibility. 
Finally,  the  constitutive  relations  for  the  viscoelastic  solid  are  given  by 


T  —  p'  I  +  2  p 
T  =  AI  +  2  p 


dip 


dIBP(t) 
dip 


+  7 


dip 


OIbo 


+  Ibg~ 


9IIBpW 
dip 


B 


dllBc 


B 


G 


'p(t) 


dip 


dip 


Oils 


B 


OIIbc 


B? 


p(() 


+ 


p(t) 


9in 


9 Df 


Dr 


Qtjm 

onG 


q,,  =  -k{0)  grad0,  s  = 

2.2.2  Specific  case 

Specifically,  we  choose  the  stored  energy  as 


*P  =  A°+(Bs+4)  (0  -  9s)—±  (0  -  0sY-40 in  (  -  )  + 


0  \  ,  PGQ  -  0G\0 

2  p0s 


(Ibg~  3)+ 


Mpo  ~  BP\0 
2  p0s 


(31) 

(32) 

(33) 

(34) 


(35) 


(IbpU,~  3), (36) 
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where  pGo- ppl  are  elastic  constants,  0S  is  a  reference  temperature  for  the  viscoelastic 
solid  and  the  rate  of  dissipation  as 


where  t 7  is  the  viscosity. 


Cm  =  r}(e)  (DG.Bp(#)DG)m, 


(37) 


Now, 


s  — 


dip 

~oo 


-  -(«*  +  d)  +  <»  -  *.)  +  din  ( f )  +  d  +  gUrVBo  -  3)  +  -  3). 


W. 


e  =  ip  +  9s 


=  a*  -  B-e,  +  <4  (#-*.)  +  a  («*  _  +  £g.(/Bo  —3)  +  -  3). 


2p9s 


hpO 


2p0s 


and  the  specific  heat  capacity  is 

Then,  (35)  reduces  to 

T  =  pi  +  pPBp(t), 


=  -e6>  =  c?0  +  c|. 


T  =  AI  +  pgBg  +  ^  (Dg-BpPjDg)"1  1  (Bp(()Dg  +  DcBp(t)) , 


n  =  W  W6  jjG  =  1‘Gi-l-aoe 

where  1 


.  From  (43) 

V 


(p  -  A)I  +  ppBp(t)  =  pgBg  +  2  (Dg-Bp^Dg)'"  (Bp(t)DG  +  DGBp(f)) , 

g-i 

and  so  by  pre  -  multiplying  the  above  equation  by  p(0  and  taking  the  trace,  we  get 


(p-A)  = 


MgMB  d,1,Bg)  —  3pp 


(38) 

(39) 

(40) 


(41) 

(42) 


(43) 


(44) 


(45) 


Using  (45)  in  (44),  we  arrive  at  the  following  evolution  equation  for  the  natural  configuration 


PgM Bp(]}Bg)  -  3pp 


I  +  PpBp(t)  =  pgBg  +  ~  (DG-Bp(()DG)m  1  (Bp(t)DG  +  DcBp(f)) 


This  is  the  evolution  equation  of  the  natural  configuration.  Thus,  with  the  current  choice  of 
energy  storage  and  rate  of  dissipation,  we  arrive  at  the  following  constitutive  equations: 

T  =  pi  +  PpBpp) 


(46) 


(47) 
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and 


HGtr(Bp{t)BG)  —  3/if 


I  +  MpBp(t)  =  ^gBg  +  ^  (Dc-B^Dg)’"  1  (Bp(t)DG  +  DGBp(tj) 


along  with 


1-r-l  V 


T~  1 


Dg  =  -2Vs(t)  Bp(t)  VKp(() 


(48) 


(49) 


The  above  model  reduces  to  the  Maxwell  fluid-  like  model  derived  by  Rajagopal  and  Srinivasa 
(2000)  when  m  =  1  and  ftc,  =  0. 

2.2.3  Application  of  Model 

We  shall  apply  this  model  to  a  deformation  that  is  homogeneous  and  in  uniaxial 
compression/extension,  given  by 


x  =  \(t),  y  = 


z  = 


y/m'  ^W) ■ 


The  velocity  gradient  is  given  by 


L  =  diagU’-4"^  • 


(50) 


(51) 


We  shall  assume 


So, 


Bp(t)  —  {B. 


1  1 


Vb'  \Jb\  ' 


B 


Bp(0  -  diaS  j5’  2Dt/2  ’  253/2  }  ’ 


v  I  •  2BX  B  X 

Bp(t)=  diag  <  B - — , + 


B  A 
+ 


A  ’  2B:i/2  X^/b'  2 £3/2  Av/B 


(52) 

(53) 

(54) 


and 


vKp(  =diag(^,  1 


1  (  'Si/4’51/4 


(55) 


_  1  ,  If?  2A  A  f?  A  5 

Dg  =  —  2diag \b-T'X-2B'X~2B 


(56) 


Also, 


g  =  v-(  f 

Kp(t ) 


=  diag  • 


A  fl1/4  f?1/4 

JX 


(57) 
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which  gives 


J  A2  \[B  \[B 

Bo  =  d,ag  B’  —  T 


and 


.  „  f  A 2  D  D 

B-,1Bc  =  d,a6|^,T.x 

Substituting  (52),  (54),  (59)  in  (46)  and  for  m  —  1 


B  B\  flG  A2 

ri+?«  „ 


AC  (1  +  2 G3/2)  J  ' 

In  general,  we  get  the  following  evolution  equation  for  the  natural  configuration 


A  -  A 


B 
2 B 


A2 


m—1  1  FG'  ft  VpB 


'pG  (A3  +  253)  --3fipB2X 
B\(1  +  2B3/2) 


Now,  plugging  (52)  into  (47)  and  using  the  fact  that  lateral  surfaces  are  traction  free 


Tn  ~  Flp  (B  s/b)  ' 


(58) 


(59) 


(60) 


(61) 


(62) 


2.2.4  Comparison  with  experimental  creep  data 

For  the  loading  process,  with  known  constant  applied  stress  and  material  properties,  (62)  was 
first  solved  for  B(t).  Then,  (60)  was  solved  with  initial  condition  ^(0)  —  fi(O)  For  the 
unloading  process,  ^  ' 1  was  set  to  zero  and  ^  was  evaluated  using  (62).  Then  using  the  initial 
condition  )  =  )  ~  \  where  is  the  time  when  unloading  is  started,  ^  is  the  instantaneous 

elastic  stretch,  ^ 1 )  during  the  unloading  process  is  evaluated  using  (60).  All  the  calculations 
were  done  in  Matlab  using  the  “ode45”  solver. 
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Figure  1.  Comparison  of  the  model  with  experimental  data  for  creep  (Falcone  and 
Ruggles  -  Wrenn  (2009))  of  PMR  -  15  at  288  °c. 

The  parameter  values  used  were  f'G  =  4  x  10s  Pa^  p,p  =  4  x  10”  Pa s  //  =  6.436  x  w'-  Pa ,  m  =  l 


The  creep  data  from  our  model  was  compared  to  the  experimental  creep  data  of  Falcone  and 
Ruggles  -  Wrenn  (2009)  for  PMR  -  15  at  288°C.  ‘fminsearch’  function  (which  uses 
Nedler  -  Mead  simplex  method)  was  used  to  minimize  the  error  defined  by 


error  =  w  x 


'E  iptheo, 


load 


- exp 


XoadY 


y  'j  exp, load ) 


+  (1  -  w)  X 


(^theo, unload  ^ exp, unload) 

E  (£earp  ,  unload ) 


(63) 


where  ethco  denotes  the  theoretical  strain  values,  eexp  denotes  the  experimental  strain  values,  the 
suffixes  load,  unload  denote  the  values  during  loading  and  unloading  processes  respectively,  w 
is  weight.  A  weight  of  w  =  0.75  was  used  since  there  are  fewer  data  points  for  the  unloading 
process.  The  best  -  fit  values  of  the  parameters  were  found  to  be  p,G  =  fip  =  4  x  108  Pa, 

77  =  6.436  x  1012  Pa. 
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3.0  DIFFUSION  OF  A  NEWTONIAN  FLUID  THROUGH  A  VISCOELASTIC  SOLID 


In  this  section,  we  shall  extend  the  model  developed  in  the  previous  section  to  include  diffusion 
of  a  Newtonian  fluid  (like  moisture). 

3.1  Preliminaries 


Let  us  consider  a  mixture  of  a  Newtonian  fluid  and  a  viscoelastic  solid.  We  shall  also  assume 
co-occupancy  of  the  constituents  i.e.,  at  each  point  x  in  the  mixture  at  some  time  t,  the  two 
constituents  exist  together  in  a  homogenized  fashion  and  are  capable  of  moving  relative  to  each 
other.  We  shall  denote  the  quantities  associated  with  the  compressible  fluid  through  the 
superscript  /  and  with  the  superscript  s  for  that  of  the  solid.  Now,  we  shall  define  the  motion  x‘ 
for  the  i  -  th  constituent  of  the  mixture  through 

x  =  x!(X',t),  i  =  f.s  (64) 

where  X'  is  the  material  point  of  the  i  -  th  constituent  in  its  reference  configuration.  We  shall 
assume  that  the  mapping  xl  is  sufficiently  smooth  and  invertible  at  each  time  t.  The  velocity 
associated  with  the  i  -  th  constituent  is  defined  as 


and  the  deformation  gradient  through 


(65) 


(66) 


The  density  p  and  the  average  velocity  (also  known  as  barycentric  velocity)  v  of  the  mixture 
are  defined  by 

P  =  v  =  ~Xyv'-  (67) 

i  i 


We  define  the  following  derivatives  for  any  scalar  quantity  </>*  by 


II 

£1* 

dl<t>1  _  d(f>‘  ■  _  d<f)‘ 

dt  dt  Ox 

dcf)1  d(f>1 

oxi  ~  ox1' 

(68) 

where 

(69) 

Hence, 

d‘f  =  9f  +  V0'.v‘, 
dt.  dt 

(70) 

and  we  shall  also  define  the  following 

r\  1 

Sflt: 

II 

05 1  05 

+ 

<1 

■«- 

< 

(71) 
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The  velocity  gradient  for  the  i  -  th  component  L'  and  the  velocity  gradient  for  the  total  mixture 
L  are  defined  by 


L'  =  Vv\  L  =  Vv. 


+  pV  +  p'b'.v'  +  E‘  +  m'.v\ 


(72) 

(73) 


where  fl,  C1  ,  r  are  the  specific  internal  energy,  heat  flux,  radiant  heating  associated  with  the  * 
-  th  component  and  El  is  the  energy  supplied  to  the  *  -  th  constituent  from  the  other 
constituents. 


Now,  taking  the  scalar  multiplication  of  Eqn.  (78)  and  v*  and  subtracting  the  resulting  equation 
from  Eqn.  (82),  we  arrive  at 


dd  .  .  ... 

pl  —  =  T'.L'  -  divq'  +  p'r‘  +  E\ 
at 


(74) 


Now,  using  f‘  ~  tP'  +  Oif ,  where  d'  are  the  Helmholtz  potential  and  specific  entropy  of  the 

-  th  constituent,  with  ^  being  the  common  temperature  of  the  constituents  at  a  point  in  the 
mixture,  Eqn.  (83)  along  with  Eqn.  (76)  results  in 

|  (pW)  +  div  (pV v‘)  =  jT’.L1  -  div  (£)  -  iq'.VO  +  ipV  +  \e‘ 


o  j  o 2 


p>  (dip1  ,d0 

J  \  ~dT  + 11  ~dt 


i  _  dj>1 

Now,  using  the  fact  that  '  <w  ,  we  can  establish  the  following  result: 

dlibl  ;  d‘ 0  d‘ip *  dip1  dl0  ( Op)1  dip1 09 \  /_  ■  dip1 

dt  1  dt  dt  00  dt  \dt  00  Ot.J  V  00 

dip1 
dt 


6=constant 


We  shall  define 


Using  the  relation  (85)  in  Eqn.  (84)  and  summing  over  *,  we  get 


I(x>')+  (£  ^Vv)  =  i  E  T-L'  -  div  (3)  -  i,W  +  ,  Q 


i  1  j  ( 


dt 


0=constant 

Equation  (87)  is  the  balance  of  entropy  with  the  rate  of  entropy  production  C  being 


*  lypi(^L 

O^1  V  dt. 


0=constant 


(75) 


(76) 


(77) 


(78) 


(79) 
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We  shall  assume  that  the  total  entropy  production  can  be  additively  split  into  entropy  production 
due  to  thermal  effects  i.e.,  conduction  (Cc),  and  entropy  production  due  to  internal  working  and 
mixing  ( Cm).  We  shall  require  due  to  the  second  law  of  thermodynamics  that  each  of  these 
quantities  be  non  -  negative.  This  implies  that 


Cc 


q.V0 


>  0, 


(80) 


^E^E^E^) 


>  0. 


(81) 


Q=constant 


We  shall  choose  q  =  —k(p,  0)V0,  k  >  0,  so  that  (89)  automatically  satisfies.  Also,  if  we  define 
the  rate  of  dissipation  :=  0(m,  then 

T'.L  ‘  +  Y,El-Y,Pi(t-^P\  =  C  >  0.  (82) 

;  ;  V  *  /  9=constant 


Assuming 


J2  E‘  +  m'v*  =  °> 


(83) 


Eqn.  (91)  can  be  re  -  written  as 


ET '■  Emv 


d1^ 

~dT 


$=constant 


=  Cm 


Now, 


=  |  (e^J  +div  (E^vj 

=  pc^ + div  (e  p'*'  (v'  - v)  )  • 


where  ip  :=  4  is  the  average  Helmoltz  potential  of  the  mixture. 

Finally,  from  eqn  (94)  and  Eqn.  (93),  we  arrive  at 

ET'Li-Em'V- 


W  ,  y 

pm+a"' 


(v'-v) 


=  Cm. 


J  0=constant 


(84) 


(85) 


(86) 


Assuming  that  all  the  components  have  the  same  Helmholtz  potential  (95)  reduces  to 


-  Cm 

9— constant 


(87) 


3.2  Constitutive  assumptions 

Let  nt  denote  the  current  configuration  of  the  mixture  and  let  kr,  nr  denote  the  reference 
configurations  of  the  solid  and  the  fluid  respectively.  Also,  let  Kp^  denote  the  natural 
configuration  of  the  viscoelastic  solid  (see  Figure  2).  Now,  if  F\  i  =  f,  s  is  the  gradient  of  the 
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reference  configuration  of  the  solid 


current  configuration  of  the  mixture 
Kt 


Figure  2.  Constitutive  Assumptions 

vs  fvi  A  F's 

motion  (usually  known  as  deformation  gradient)  '  v  ’  /  and  Kp(‘>  is  the  gradient  of  the 
motion  of  the  viscoelastic  solid  from  kp(0  to  Kt,  then 


G"  :=  Fs 


(88) 


We  shall  assume  that  the  specific  Helmholtz  potential  for  the  mixture  is  of  the  form 

^  ^  i6'  /bo’//b«’///bg’  lBU.t)'IlBUt)'IIlBh*))  ' 


And 


<H’  ds’ip 

-  =  — +  (v-v).grad^ 


clip 

dt 


dip  dip 


dip 


dIIB 


B 


P(i) 


P(<) 


B 


p(t) 


(89) 


(90) 

(91) 
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+ 


dip  dip 

+  JSCQ 


QIBg  dll. 


dip 


Ba 


dIIBc  °  dIIIBn 


B 


-l 


Be 


.bg  +  +  {v  -  vs 


where  we  have  dropped  the  superscript  s  over  B<?'  Bp(0  for  convenience,  and  (  )  is 
Now, 


F  =  FKp(t)G  +  FKp(t)G 


=>■  FF_1  =  Fk  ,  ,GG-1F“1  +  Fh.  .  G 

rip(<)  Kp(o  p(o 

=»  L  =  Lp(t)  +  FKp(t)LGFK 


p(0 


D  -  Dp(d  +  2  (F«p(0 LgFS(0  +  fkJwLgfkpW)  > 


where  (•)  is  the  time  derivative  of  the  second  order  tensor.  In  addition, 

BpW  =  F,p(0Fr  +  FF^p(() 

=  Lp(t)Bp(t)  +  B/»(t)Lp(t), 

and  similarly 

Bg  =  LgBg  +  BgLg. 

Hence,  from  (102)  and  (103),  we  have 

Bp(0  =  LBp(0  +  BP(0L,T(t)  -  F/tp(f)  (LC  +  L£)  F^t), 

Based  on  the  assumption  of  isotropic  elastic  response,  we  choose  kp( 0  such  that 

F  =  V 

Kp(t)  Kp(0’ 

V  F 

where  Kp(‘>  is  the  right  stretch  tensor  in  the  polar  decomposition  of 

Using  (102),  (103),  (104),  (106)  in  (101),  we  get 


+2 


ip  =  2 
dip 


dip  |  r  dip 

dhZ+  B««dlu 


B 


'p(t) 


dip 


B 


dip 


dlln  “lp0) 

^pd)  vinvp{l) 


v (0  , 


(Ls 


+  / 


dl Bq  ’  ~BGdIIBG 

Volume  additivity  constraint: 


Bg 


G  ~  ~ 


dj>  2  _ 
dIIBa  G  dlllBa 

(ps  -\-  (pf  —  1, 


CM  ‘ 

•LG  +  50  ^  +  (v  ~  v' 


).grad,0, 

(92) 

dpn 

dt  . 


(93) 


(94) 


(95) 


(96) 


(97) 


-L£) 

(98) 

*).gradV’, 

(99) 


(100) 
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where  Pr,  *  —  s-  /.  Also,  we  shall  assume  that  the  constituents  are  incompressible  in  their 
natural  states,  i.e.. 


plR  =  constant. 

Eqn.  (76)  using  (110)  can  be  re  -  written  as 

<9  <j>1 


(101) 


Ot 


+  div  =  0, 


(102) 


which  implies 


^+div(x:^vj=o 

=>  div  ^e//v-^  +  </>s- vs^  =  0(using  (109)). 

Eqn  (113)  can  be  re  -  written  as 

</>str(Ls)  +  <//  tr(L-l)  +  vhgrade/f  +  v^.grade//  =  0 
Again  from  (109),  we  have 

grade/)3  +  grade//  =  0, 

Hence,  using  (1 15)  in  (1 14),  we  arrive  at 

(j>atr(L8)  +  (j)ftr(Lf)  +  vsj. grade/)3  =  0, 

where  v-s-/  =  v  —  v/,  is  the  velocity  with  which  the  fluid  diffuses  with  respect  to  the  solid. 
We  shall  also  assume  that  the  rate  of  entropy  production  is  of  the  fonn 

Cm  =  Cm  (l<G)  Fp((),  i/,  0,  Vsj)  . 

Now  (96)  along  with  (80)  reduces  to 


(103) 

(104) 

(105) 

(106) 

(107) 


(108) 


T3.L3  +  Tf.Lf  -  m3o 


S ,/  P 


M\ 


=  U(l 


dt  )  6=constant 

Using  (108)  in  (118),  we  get 

T*  L*  +  Tf.Lf  -  -  TJ,,,.  (L*  -  LJ)  -  Tfc.LJ,  -  p(v  -  v*).  (grad(/>)#_ 

=  fm(L&.F^),L^0,v.,/) 

where 


(109) 


constant 


(HO) 


Tp(0  := 


Oxp 


01 B 


+ 1 


D 


04’ 


p(t ) 


p(t ) 


01 1 B 


pit) 


_  0 E2 

m  anB,m  "<■) 


04 ) 


dIIIBpW 


(HI) 
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Tfc  :=  2p 


*  la,- 


B  -G- 


dxp 


diiBc  L'  diiiBc 


constant 


diBc  aadiiBc 

Eqn.  (119)  with  the  constraint  (116)  can  be  written  as 

T\LS  +  Tf.V  -  m *.vsJ  -  Tsp{t).  (L6  -  Lfc)  -  Tfc.Lfc  -  p(v  -  Vs).  (grad^)*=. 

+  A  j>tr(L6)  +  <^tr(l/)  +  vS)/.grad^)  =  £m  (l/G, F6(t),  i/,  0,  vsJ) 

We  shall  further  assume  that  the  rate  of  dissipation  can  be  additively  split  into  the  rate  of 
dissipation  due  to  mechanical  working,  rate  of  dissipation  due  to  the  fluid  and  the  rate  of 
dissipation  due  to  diffusion  of  the  fluid,  with  specific  forms  as  follows: 

U  (lg.  Fp(t),l/,  0.  va>/)  =  £  (Lg,  Bp(t),  0)  +  +  a(0)vsJ.vsJ. 

Then,  from  (123)  and  (122),  we  arrive  at 

T  Ls  +  Tf.L-f  -  ms.vsj  -  T‘,t).  (L6  -  LG)  -  T°G.LSG  -  p(v  -  v6).  (grad^^ 


(112) 


(113) 


(114) 


stant 


+ 


A  (Vtr(L6)  +  cfS  tr(Lf)  +  vsJ.grad</>s)  =  £  (LG,  Bp(t), 6)  +  v. D/.D/  +  a(0) 


W,/-vs.f 


(115) 


which  can  be  re-  written  as 

L”.  [t*  +  M‘l  -  TJ(t|]  +  \J .  [t'  +  M/1  -  •'V'}  +  (t;w  -  T-o)  .Li 


G 


+  vsJ.  [-ill"  +  A  graded  -  a(9)vsJ  +  pf  (gradt/>)*=conatant  =  £  (Lsa.Bp{t),6)  , 

using  the  fact  that  P  fv  —  v  )  =  Pfvsj.  Since,  the  right  hand  side  of  (125)  does  not  depend  on  L 
,  B- 1  and  V«J,  we  have 

Ts  =  -A0SI  +  T*(t), 

Tf  =  +  uDf, 

ms  =  A  grade))6'  -  a{6)vsJ  +  pf  {grad^)8=constant , 

and  so  (125) reduces  to 

(T;W-T6G)  .L6G  =  £(L^B,,(t).0) 


(117) 

(118) 
(119) 


(120) 

(121) 


Now,  we  have 

pr  =  psjs  =  psJgj;, 

where  ^  =  det(Fk),  k  —  G.p  or  no  subscript.  From  (110)  and  (130)  and  further  assumption 
that  there  is  no  diffusion  of  mass  from  kp(0  to  Kt,  we  have 

dsJr 

=  o  =>  trLfc  =  0. 

c  It  G  (122) 


Jq  =  constant 
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Next,  we  shall  maximize  the  rate  of  dissipation  £  with  (129)  and  (130)  as  constraints.  We  shall 
maximize  the  auxiliary  function 


Now, 


*  :=  ■£  +  0  [t  -  (T;(f)  -  Tb)  -L h\  +  H  (I.L aG) 

_  (p  + 1)  dt  /  \  Mt_o 

r  -  (,TP(0  -  T a)  +  01-0- 


OLb 


=  0 


P  dL 


G 


Taking  the  scalar  product  of  (133)  with  lg  and  using  (129),  (131) 

(0  +  1)  £ 


0  JgW 


Next,  taking  the  trace  of  (133) 

/f  =  1 

0  3 


The  evolution  equation  for  the  natural  configuration  of  the  solid  is  given  by 


(^)-T&)-|l+1 


(123) 


(124) 


(125) 


(126) 


(127) 


We  shall  assume  the  following  specific  form  for  the  specific  Helmholtz  potential  of  the  mixture 


4>  =  A'Mir+4)  (« -  «,)-f  (#  -  ».?-<%>  in  ({)+— <.,~3)- 

(128) 

where  ^Go?MGi)/h>o>  MPl  are  elastic  constants,  ^  is  a  reference  temperature  for  the  viscoelastic 
solid 

Now, 

s  =  -00 

= -(B* + « + ■*  «*  -  ">> + ^in  (f  H + - 3) + -*  (129) 

The  internal  energy  e  is  given  by 

e  =  ip  +  9s 


=  A‘-  art.  +  4(t> - «.)  +  |  (#2  - #;)  +  ££</*,  - 3)  +  0-(/flp(el  - 3). 


and  the  specific  heat  capacity  C7  is 


From  (137)  and  (120) 


Cv  =  -eO  =  cl9  +  4. 


(130) 

(131) 
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(132) 


rps  P  W& 

p{t)  P°R  Bs 

and  from  (137)  and  (121), 


■  -0f-  K-  - 3) 1  -  ~ 3>  >]  • 


Ts  —  ^  R  _  Pp  ( j  _  o\  x 

G  psR  9S  ['  G  G  2(J£)2  YD^  ) 


1?T2  Uflo  -  3) 1 


2  (Jsg) 


where  Pg  —  Pgo  ~  Pgi@,  Pp  ~  PpO  Ppi@. 

We  shall  further  assume  that  the  rate  of  dissipation  £  is  of  the  form 

£  =  7(0)Dg.Dg, 


then  (126)  becomes 


T^-A^  +  ^-r  fipBp(t) 


fopm  -  3) 


^2^G-3)I 


2  (Ji) 


and  (136) reduces  to 


| PpBp(t)  ~  Pg'Bg  ~  £  [/h>tr(Bp(f))  -  Mctr(BG)]  1 1  =  7(0)Dg 
The  final  constitutive  equations  are 

T'  =  -w  +  £  [ft>Bp«)  -  T.  -  3)  1  -  17  (’Bo  -  3)  I 

r  L  ,JP  X  7  up 

Tf  =  -\4>fI  +  vBf, 

ms  =  A  graded  -  a(0)vaj  +  pf  (grad ip) e=constant , 


_  llp0  -  Oflp  1  _  _  //G0  -  Ougi 


where 


Pg  = 


0  s  with 


P-  j/JpBpW  “  Pg&g  -  ^  [pPtv(BpW)  -  /iGtr(BG)]  1 1  =  7(fl)DG 


being  the  evolution  equation  of  the  natural  configuration. 


(133) 

(134) 

(135) 

(136) 

(137) 

(138) 
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